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We show that quantum subdynamics of an open quantum system can always be described by a Heimitian 
map, irrespective of the form of the initial total system state. Since the theory of quantum error correction 
was developed based on the assumption of completely positive (CP) maps, we present a generalized theory 
of linear quantum error correction, which applies to any linear map describing the open system evolution. In 
the physically relevant setting of Hermitian maps, we show that the CP-map based version of quantum error 
correction theory applies without modifications. However, we show that a more general scenario is also possible, 
where the recovery map is Hermitian but not CP. Since non-CP maps have non-positive matrices in their range, 
we provide a geometric characterization of the positivity domain of general linear maps. In particular, we show 
that this domain is convex, and that this implies a simple algorithm for finding its boundary. 
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I. INTRODUCTION 



The problem of the formulation and characterization of the 
dynamics of quantum open systems has a long and extensive 
history iH |3, l3|]. This problem has become particularly rel- 
evant in the context of quantum information processing f?], 
where a remarkable theory of quantum error correction (QEC) 
was developed in recent years to address the problem of how 
to process quantum information in the presence of decoher- 
ence and imperfect control jH]. A key assumption common to 
many previous QEC studies is that the evolution of the quan- 
tum information processor can be described by a succession 
of completely positive (CP) maps interrupted by unitary 
gates or measurements ff]. However, it is well known that 
if the initial total system state is entangled, quantum dynam- 
ics is not described by a CP mapjl |1 [M 111 [H- In fact, 
we showed very recently in Ref. fiy\ that a CP map arises if 
and only if the initial total system state has vanishing quantum 
discord lO . i.e., is purely classically correlated. One is thus 
naturally led to ask whether this impacts the applicability of 
QEC theory under circumstances where non-classical initial 
state correlations play a role. Here "initial state" does not re- 
fer exclusively to the "t — 0" point, but also to intermediate 
times where the recovery map is applied, since this map was 
also assumed to be CP in standard quantum error correction 
theory [7J. Motivated by this fact we here critically revisit the 
CP maps assumption in QEC, and show that it can be relaxed 
' . To do so, we first consider the problem of characterizing 
the type of map that describes open system evolution given 
an arbitrary initial total system state (Section HH). We show 
that this map is always a linear, Hermitian map (of which 
CP maps are a special case). We then argue that the generic 
noise map describing the evolution of a quantum computer 
as it undergoes fault tolerant quantum error correction (FT- 



' Note that this is issue is entirely distinct from the critique of Marko- 
vian fault tolerant QEC expressed in fT5], which was concerned with the 
compatibility of other assumptions of fault-tolerant QEC (specifically, fast 
gates and pure ancillas) with rigorous derivations of the Markovian limit. 



QEC) is indeed not a CP map, but rather such a Hermitian, 
linear map (Sectionllllb. The reason is, essentially, that imper- 
fect error correction results in residual non-classical coiTela- 
tions between the system and the bath, as the next QEC cycle 
is applied. To deal with this, we develop a generalized the- 
ory of QEC which we call "linear quantum error correction" 
(LQEC), which applies to arbitrary linear maps on the system 
(Section IIVI ). Then we show that, fortunately, the CP-map 
based version of QEC theory applies without modifications 
in the physically relevant setting of Hermitian maps. How- 
ever, we show that a more general scenario is also possible, 
where the recovery map is Hermitian but not CP. This is useful 
since it obviates the unrealistic assumption that the recovery 
ancillas enter the QEC cycle as classically correlated with the 
other system qubits. Our results significantly extend the realm 
of applicability of QEC, in particular to arbitrarily correlated 
system-environment states. We conclude in SectionlV] 



II. QUANTUM DYNAMICAL PROCESSES AND MAPS 

In this section we prove a basic new result, that a quantum 
dynamical process can always be represented as a linear, Her- 
mitian map from the initial to the final system-only state. In 
doing so we rely heavily on our previous work il3ll . 

The dynamics of open quantum systems can be described 
as follows. Consider a quantum system S coupled to another 
system B, with respective Hilbert spaces Tis and T-Lb, such 
that together they form one isolated system, described by the 
joint initial state (density matrix) pss(O). Their joint time- 
evolved state is then 



PSB 



{t)^uit)psBmHt), 



(1) 



where U (t) is the unitary propagator of the joint system-bath 
dynamics from the initial time t = to the final time t, i.e., 
the solution to the Schrodinger equation U = —{i/h)[H, U], 
where H is the joint system-bath Hamiltonian. The object of 
interest is the system S, whose state at all times t is governed 
according to the standard quantum-mechanical prescription 
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by the following quantum dynamical process (QDP): 

ps{t) = TrsbssW] = TrB[C/ss(i)pss(0)C/ssW^]- (2) 

Ttb represents the partial trace operation, corresponding to 
an averaging over the bath degrees of freedom f3^ . 

The QDP (|2]i is a transformation from psb(O) to ps{t)- 
However, since we are not interested in the state of the bath, it 
is natural to ask: 

Under which conditions on psb(O) is the 
QDP a map 'I>Q(t), 

ps{t) ^ ^Q{t)[psm, (3) 

and what are the properties of this map? 

In general, a map is an association of elements in the range 
with elements in the domain. Here we use the term "map" 
solely to indicate a state-independent transformation between 
two copies of the same Hilbert space, in particular Tis ^ 
Tis^ Then, a well-known partial answer is that if psb(O) is a 
tensor product state, i.e., psb{0) = Psi^) ® Pb{0), then the 
QDP (|2]i is a CP map. A more general answer was provided 
in [13]. To explain this answer we must first introduce some 
terminology. 

A. Various linear maps 

A map $ : B{Ti.) ^-^ B{'H) [space of bounded operators 
on Ti] is linear if $[api + bp2\ = + fc'S[/02] for any 

pair of states pi, p2 : H i-^ H, and constants a,5 e C. A 
linear map is called Hermitian if it maps all Hermitian opera- 
tors in its domain to Hermitian operators. We first present an 
operator sum representation for arbitrary and Hermitian linear 
maps, that generalizes the standard Kraus representation for 
CP maps [j6t]. The proof is presented in Appendix lAl 

Theorem 1 A map $l : SO^n ^ 9Km (where 3Jl„ is the space 
of n X n matrices) is linear iff it can be represented as 

<1>l(p) = Y^ E^pE2 (4) 

where the "left and right operation elements" {Ea} and 
{E'^} are, respectively, m x n and n x m matrices. 
$H is a Hermitian map iff 

*H(p)-y c^EapEl c„gR. (5) 

We will sometimes denote a linear map by listing its ele- 
ments, as in $L = {^a, E'^y^^^. Note that a linear map 
$L = {-Eq,£'^}J;^i is trace preserving if YJa^i^aEc = 



^ This is meant to exclude claims that system state-dependent transforma- 
tions qualify as CP maps, as in Ref. ||161 . In such cases the elements of the 
transformation (the "Kraus operators") depend on the system input state, 
which contradicts our notion of a map. 



/. Also note that the two sets of operation elements 
{E^,E'^}\^^jmd {Fp^F'^Yp^^, where Fp = ELi ^«/3^^« 
and F'p ~ Yl^a=i ^apE'^, represent the same linear map <i>L if 
the matrices u and v satisfy uw^ — /. 

As a simple example of a non-CP, Hermitian map, con- 
sider the inverse-phase-flip map. The well-known CP phase- 
flip map is 101: ^ppip) = (1 ^ p)p + P<7zP<^z, where 

< p < 1 and ffz is a Pauli matrix. Solving for $pp from 
*pf[^'pf(p)] = p, we find that $pp(p) = cop + cicr^pcr^, 
where ci — p/ {2p — 1) and cq = 1 — ci, and cq, ci have op- 
posite sign for < p < 1. Moreover, Tr[<i>pp(p)] = Tr(p). 
Therefore <i>pp is a trace-preserving, Hermitian, non-CP map. 

A linear map is called "completely positive" (CP) if it is 
a Hermitian map with Cq > Va. CP maps play a key 
role in quantum information and quantum error correction 101, 
though they have a much earlier origin There are other 

useful characterizations of CP maps - see, e.g., Refs. HH. It 
turns out that there is a tight connection between CP and Her- 
mitian maps iflolfTlll : a map is Hermitian iff it can be written 
as the difference of two CP maps. 

The definition of a CP map $cp implies that it can be ex- 
pressed in the Kraus operator sum representation |6i]: 

ps{t) - ^i?„(i)p5(0)i?i(i) = <f>cp(t)[ps(0)]. (6) 

a 

If the operation elements satisfy ^^El^Ea = I then 
Tr[p5W] - 1. 

B. Special linear states 

Following Ref. fisl], we define the class of "special-linear" 
(SL) states for which the QDP ^ always results in a linear, 
Hermitian map. An arbitrary bipartite state on Tis <E)'Hb can 
be written as 

P5S = V . . Qij \i) {j I «) (jjtj , (7) 

where is an orthonormal basis for Tis, and 

{4>ij}i%'i'^ '■ T^B ^ ai-e normahzed such that if 

Tr[0ij] ^ then Tr[(/<ij] = 1. The corresponding reduced 
system and bath states are then ps = j)£C 
where C = {(i, j)|Tr[0ij] = 1}, and pb(0) = Xli £'"<^"- 
Hermiticity and normalization of psB, Ps^ and pB imply 
eij = Pji' '/'y = and J2i Qti = 1- 

Definition 1 A bipartite state psB, parametrized as in Eq. 
0, is in the SL-class if either Tr[(/)jj] ~ \ or (pij — 0, Vi, j. 

Thus a non-SL state is a state for which there exist indexes 

1 and j such that T^r[4>ij]_^= but (j>ij ^ 0. The following 
result proven in Ref. [13] (generalizing an earlier result in 
Ref. UM) provides an almost complete answer to the question 
posed above: 

Theorem 2 (Theorem 2 of illj) // P5_b(0) is an SL-class 
state then the QDP (|2| is a linear, Hermitian map <I>h : 
PsiO)^Psit). 
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A further result proven in Ref. ifTsIl (Theorem 3 there) pro- 
vides necessary and sufficient conditions on Psb{0) for the 
QDP (|2]i to be a CP map, namely, p<; b (0) should be a state 
with vanishing quantum discord [14]. Such a state cannot 
contain any quantum correlations. This clearly illustrates the 
limitations of CP maps in describing quantum dynamics. At 
the same time one may wonder as to the generality of the 
SL-class employed in Theorem |2l Non-SL states are sparse 
ifisll . so it is in this regard that we stated that Theorem|2pro- 
vides an almost complete answer to the question posed above. 
However, we can go further. As mentioned without proof in 
Ref. 1 13], in fact the QDP (|2]i is a linear, Hermitian map from 
Ps(0) Ps{t) for any initial state pss(O). We next prove 
this key fact. 

C. Hermitian maps for arbitrary initial states 

We split the general initial state representation O into a 
sum over SL and non-SL terms (thus splitting { Qij } and {0^^ } 
into two sets): 

PSb{0)= Y1 (^^]\^){j\®V^3 + f^vl'^) (jl ^ ■ 

iiG(SL) iie(nSL) 

(8) 

In accordance with the definition of SL states, in the first sum 
we include only terms a,;j |i) {j\ (pij for which Tr[(p,;j] ^ 
or (fiij = 0, in the second only terms {j\ ® i/jij with bath 
operators satisfying ^ and Tr[-i/;ij] = 0. By virtue 
of this decomposition only the first term contributes to the ini- 
tial system state: ps{0) = Trs[psB(0)] = J2^j{sh) "»jl*>0'l- 
This is because the condition Tr[-0y] — eliminates any con- 
tribution from the second term in the decomposition (O to the 
initial system state. Consequently Eq. Q assumes an affine 
form: 

*qW[P5(0)] = $SL(i)[ps(0)] + KnSL(i), (9) 

with the term KnShit) being a shift that is independent of 
Ps{0). 

As shown in Ref. 01311 . the linear map $gL is constructed 
as a function of the bath operators {fij}'. 

^sLmpsm^ J2 x:iv,%p.psio)p,{w^,,)\ 

(ij)e(SL);fc,a 

(10) 

where Pi = are projectors, X^J are the singular values 
in the singular value decomposition 0^ = J^a X^i \xfj){yij l> 
and the operators V^^ = {^Jk\U\xf^) and W^^^ = {AlUlyfj) 
act on the system only, with {[V^fc)} being an orthonormal ba- 
sis for the bath Hilbert space Hb- 

In addition, the non-SL terms in Eq.® generate the shift 
term 

KnSLit)^ P^J^TB[UsBm){J\(E>^P^JUlgit)]. 

(11) 

This shows explicitly that KnSh{t) does not depend on the 
initial system state, since the latter is fully parametrized by the 



coefficients {ctij}ij £{sl)^ while K^shit) depends only upon 
the coefficients {Aj}ijg(nSL)- 

Now we take a further step to argue that the affine map (|9]l 
is actually a linear, Hermitian map if the map acts only on the 
space of density matrices. This is a direct application of the 
result in Ref. iflOl . 

Theorem 3 The QDP (|2) is representable as a linear, Hemii- 
tian map ^ii{t) : ps{0) i-^ ps(t) for any initial system-bath 
state. 

Proof. Let N = dimHs- Let Fq = I and let {F^ : Tr(F^) = 

0}^^i^ be a basis for the set of traceless Hermitian matrices 
which are mutually orthogonal with respect to the Hilbert- 
Schmidt inner product, i.e., Ti{F^F^) = NSf^.^,. Hence the 
initial system state ps{0) can be expanded as 

ps{0) - -(/+ J2 = Tr[ps(0)F^] = (F^)p,(o), 

(12) 

and the final system state is found to be 

Ps{t) = ^[<i>SL(/)+ E &M*Sl(^J] + i^nSL 

= ^Yi{t)[psm, (13) 

where the equivalent Hermitian map is constructed by set- 
ting $h(/) = *sl(/) + NK^si. and $h(Fm) = ^SL(i^;.) 
1 < V/i < N'^ — 1. That this map is Hermitian is simple to 
verify, for all the components are Hermitian. ■ 

TheoremOprovides a complete, and perhaps surprising an- 
swer to the question posed at the beginning of this section. 
Namely, the most general form of a quantum dynamical pro- 
cess, irrespective of the initial system-bath state (in particular 
arbitrarily entangled initial states are possible) is always re- 
ducible to a Hermitian map from the initial system to the final 
system state. The surprising aspect of this result is that it was 
not known previously whether QDP could always even be re- 
duced to a map between system states. 

Of course, this result does not resolve the more difficult 
question of ensuring the positivity of the final system state. 
That is, a Hermitian map may transform an initially positive 
system state to a non-positive one, violating the postulate of 
positivity of quantum states. To resolve this one must identify 
the "positivity domain" of $h, i-C-, the set of initial system 
states (positive by definition) which are mapped to positive 
states by $h ifToll - We address this in the next subsection. 

D. Geometric characterization of the Positivity Domain 

In this subsection we prove the convexity of the positivity 
domain and propose a geometric method for characterizing it. 
Let S{n) = {/9 e CiJ-L) : p > 0,Trp = 1}, where C{n) is 
the set of all linear operators on H. The positivity domain of 
a linear map $l : S{n) i-> B{n) is: P^ = {p e 5(Hs) : 
<i>L{p) > 0}. 
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Following earlier work iH [Tl, H^l, in Ref. (tIS, a com- 
plete geometric characterization of density matrices was given 
by using the Bloch vector representation for an arbitrary A^- 
dimensional Hilbert space H. This works as follows: let 

{F^}'^=i^ be a basis set as in the proof of Theorem |3] 
whence the expansion (fTZb applies again. The vector b = 
(61, ...,6jv2^i) S M.-^ ~^ of expectation values is known as 
the Bloch vector, and knowing its components is equivalent 
to complete knowledge of the corresponding density matrix, 

via the map b 1-^ p=jj{I + 'Y^^=i^ b^Ffj_). Let n denote a 
unit vector, i.e., n G and X^iLi^^^f — 1' de- 

fine F„ = X]^=r^ "-M^M- L^*- '■^^ minimum eigenvalue of 
each F„ be denoted m(i^„). The "Bloch space" B{R^^-^) 
is the set of all Bloch vectors and is a closed convex set, since 
the set S{H) is closed and convex, and the map b 1— > p i s 
linear homeomorphic. As shown in Theorem 1 of Ref. luHl . 
the Bloch space is characterized in the "spherical coordinates" 
determined by {F„} as: 

B(R^'-i) = (b = m£R^'-i:r< , } , 1. (14) 
I \m{Fn)\} 

It is hard to imagine a more intuitive or simpler geometric 
picture. 

Next we show that the positivity domain is a convex set as 
well. 

Proposition 1 The positivity domain P$ of a linear map $l 
is a convex set. 

Proof. Consider two density matrices p and p' as inte- 
rior points of P$ with corresponding Bloch vectors b — 
(61, 6Ar2_i) and b' = (6']^, 5^2_i)- The claim is that 
a third density matrix p" with corresponding Bloch vector 
h"{a) = Q!b+ (1 — Q;)b', with <a< 1, is then also interior 
to P$. This follows directly by linearity of the map <i>L. First, 

by assumption $l[/o] = + b^F^,)] > and 

<fL[p'] = + W] > 0' so that + 

(1 - a)$L[p'] > 0. Second, a$L[p] + (1 - a)<Pi^[p'] = 

+ E^=r' ^ Therefore indeed 

> 0. ■ 

We are now ready to describe an algorithm for finding the 
boundary of the positivity domain P$. We know at this point 
that P$ is convex and that P$ is a subset of the Bloch space, 
itself a closed convex set. Pick a unit vector n and draw a line 
through the origin of the Bloch space along n. If P$ includes 
the origin, i.e., the maximally mixed state, then convexity im- 
plies that this Une intersects the boundary of P$ once. If P$ 
does not include the origin then convexity implies that this 
line either intersects the boundary of P$ twice or not at all. 
I.e., it follows from convexity that the line may not re-enter 
the positivity domain once it exited. In order to determine this 
boundary we may thus compute the eigenvalues of ^Lbnl?")] 
as a function of r, where r is the parameter in Eq. ( fT4b . and 
where pn{r) is the density matrix determined via the mapping 



b = rn I— > p. The computation should start from r = and 
go up to at most r = l/|77i(P„)|. The boundary is identi- 
fied as soon as the eigenvalues of <I>L[pn(r)] go from all pos- 
itive semi-definite to at least one negative, or vice versa. For 
each unit vector n, the corresponding point on the border of 
the positivity domain can be found in this way. Then the al- 
gorithm constructs the boundary of the positivity domain by 
finding the boundary points in all directions n. Of course, in 
practice one can only sample the space of unit vectors n and 
factors r. In principle this yields a complete geometrical de- 
scription of the positivity domain of a given linear map. 



III. CP MAPS AND FAULT TOLERANT QUANTUM 
ERROR CORRECTION 

A. CP maps: pro and con 

We have already mentioned that a QDP (|2]i becomes a CP 
map iff the initial system-bath state has vanishing quantum 
discord, i.e., is purely classically correlated [13]. The stan- 
dard argument in favor of CP maps is that since the system S 
may be coupled with the bath B, the maps describing phys- 
ical processes on S should be such that all their extensions 
into higher dimensional spaces should remain positive, i.e., 
^CP (81 > Vn G Z+, where /„ is the n-dimensional 
identity operator. However, one may question whether this is 
the right criterion for describing quantum dynamics [8J. An 
alternative viewpoint is to seek a description that appUes to 
arbitrary psb{0), as we have done above. We now argue that 
this viewpoint is the correct one for fault-tolerant quantum er- 
ror correction (FT-QEC). 



B. (In)validity of the CP map model in FT-QEC 

Let us show that system-environment correlations impose a 
severe restriction on the applicability of CP maps in FT-QEC. 
The CP map model used in FT-QEC ll22ll2l [24i [Isj [2^ [iTj 
28, 29] can be described as follows (see, e.g., Eq. (8.1) in 
ii): ps{T) = $^°*(r,io)[ps(io)] where 

$^°^(r,to) = ^^u{U)^Gp{U,k^i), (15) 

l—l 

where T = is the total circuit time, and where ^(/[ps] = 
UsPsU\ is a unitary map (automatically CP) that describes 
an ideal quantum logic gate.^ This represents the idea used 
repeatedly in FT-QEC, that the noisy evolution at every time 
step can be decomposed into "pure noise" ^cp{ti, ii-i) fol- 
lowed by an instantaneous and perfect unitary gate ^u{ti). 
More precisely, in FT-QEC one assumes that the evolution 
starts (t = tQ — 0) from a product state, then undergoes a CP 



^ In this subsection we denote noise maps by their initial and final times, to 
distinguish them from the instantaneous unitary maps. 
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map ^cp 1 1 ^o) due to coupling to the environment, followed 
by an instantaneous error correction step ^i;{ti). If the latter 
were perfect then the post-error-correction state would again 
be a product state ps{ti) ^ Psiti)- However, FT-QEC allows 
for the fact that the error correction step is almost never per- 
fect, which means that there is a residual correlation between 
system and bath at ti. Hence, according to Ref. [Ts"], the map 
that describes the evolution of the system is a CP map if and 
only if the residual correlation is purely classical. Otherwise it 
is a Hermitian map. To make this point more explicit, consider 
a sequence of two noise time-steps, interrupted by one error 
correction step. In the ideal scenario, where the error correc- 
tion step ^u{ti) works perfectly (i.e., reduces the system-bath 
correlations to purely classical), we would have 

^^cU^2M = *cp(i2,ii)*t/(ii)$cp(ii,to), (16) 

where $cp(^2, ii) is again a CP noise map. However, in re- 
ality ^u{ti) works imperfectly [system-bath correlations are 
not purely classical after the action of and the actual 

map obtained is 

*L^Hi2,to) = $H(i2,ti)$c/(ti)*cp(ii,io), (17) 

where $H(i2, ti) is now a Hermitian map. Note that, in fact, 
even the assumption that the first noise map is CP will not 
be true in general, due to errors in the preparation of the ini- 
tial state, leading to non-classical correlations between system 
and bath. We conclude that in general the CP map model ( flSl l 
should be replaced by 

^'^'{TM) - (g)^^ ^u{k)^u{k,t,^{), (18) 

where $H(ii, ^i-i) are Hermitian maps, not necessarily CP.'* 
It is worth emphasizing that this distinction between purely 
classical and other correlations, and the resulting difference 
between CP and Hermitian evolution, is not a distinction that 
has thus far been made in FT-QEC theory. Rather, in FT- 
QEC one distinguishes between "good" and "bad" fault paths, 
where the former (latter) contain only a few (too many) errors. 
Quoting from |30]: "There are good fault paths with so-called 
sparse numbers of faults which keep being corrected during 
the computation and which lead to (approximately) correct an- 
swers of the computation; and there are bad fault-paths which 
contain too many faults to be corrected and imply a crash of 
the quantum computer." This leads to a splitting of the total 
map ( fTSl l into a sum over good and bad paths. One then shows 
that the computation can proceed robustly via the use of con- 
catenated codes, provided the "bad" paths are appropriately 
bounded. In 1 28] (p. 1272) it was pointed out that the sum over 
"good" paths need not be a CP map, but can be decomposed 
into a new sum over CP maps [Eq. (8.13) there]. This new de- 
composition can then be treated using standard FT-QEC tech- 
niques. However, this assumes again that the total evolution is 
a CP map, which in fact it is not [Eq. (fTsTll. 



Note that Eq. (Ts) applies also to non-Markovian noise, and is hence com- 
plementary to Hamiltonian FT-QEC 13^(31^3^. 



These observations motivate a generalized theory of QEC, 
which can handle non-CP noise maps. This is the subject of 
the next section. The main result of this theory is reassuring: 
in spite of the invalidity of the CP map model in FT-QEC, the 
CP-map based results apply because the same encoding and 
recovery that corrects a Hermitian map can be used to correct 
a closely related CP map, whose coefficients are the absolute 
values of the Hermitian map. This is formalized in Corollary 

m 

IV. LINEAR QUANTUM ERROR CORRECTION 

Having argued that non-CP Hermitian maps arise naturally 
in the study of open systems, and in particular FT-QEC, we 
now proceed to develop the theory of Linear QEC. For gen- 
erality we do this for arbitrary linear maps, i.e., maps of the 
form (|4|i. We then speciaUze to the physically relevant case of 
Hermitian maps. 

Let us first recall the fundamental theorem of "standard" 
QEC (for CP noise and CP recovery maps) ff]: Let P be a 
projection operator onto the code space. Necessary and suf- 
ficient conditions for quantum error correction of a CP map, 
<fcp(p) = E.^*P^/ai-e 

PF}FjP^\^jP (19) 

An elegant proof of this theorem and a construction of the 
corresponding CP recovery map was given in Refs. ]4, 33]; 
we use some of their methods in the proofs of Theorems 1415 1 

A. CP-recoverable linear noise maps 

While general (non-Hermitian) linear maps of the form (|4| 
do not arise from quantum dynamical processes [Eq. (|2|l], it 
is still interesting from a purely mathematical standpoint to 
consider QEC for such maps. Moreover, we easily recover 
the physical setting from these general considerations. 

Theorem |4] shows that there is a class of linear noise maps 
which are equivalent to certain non-trace-preserving CP noise 
maps when it comes to error correction using CP recovery 
maps. 

Theorem 4 Consider a general linear noise map $l (p) — 
"^i^iEipE'^ and associate to it an "expanded" CP map 
^cAp) = hT.tlE^pEl + \Y.tiE[pK^- Then any QEC 
code C and corresponding CP recovery map TZ for <1>CP ore 
also a QEC code and CP recovery map for <i>L. 

Proof. The operation elements of <&cp are {Fi}fLi = 
{^E,}l^ and {F^+..}fL, = ij^mti^ whence 
'^cp(p) — Y^^'i^iEipF^ . The standard quantum error cor- 
rection conditions (fT9b for <E>cp, where 

A-2(;t J,)=At, (20) 
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become three sets of conditions in terms of the Ei and i?': 

(i) PEjEjP = 2a,, P, (ii) PE'^E'^P = 2a[,P, 
(in) PEJe'jP = 27yF, (21) 

where i,j G {1,...,7V} and Uij = Ay, 7^ = Xi,N+j, 
ctij = XN+i,N+j- The existence of a projector P which sat- 
isfies Eqs. (|2T]i(i)-(iii) is equivalent to the existence of a QEC 
code for $cp- Assuming that a code C has been found (i.e., 
PC = C) for $cp^ we use this as a code for $l ™d show that 
the corresponding CP recovery map TZqp is also a recovery 
map for Indeed, let Gj = Y^'^=i '^ijPi be new operation 
elements for $cp^ where u is the unitary matrix that diago- 
nahzes A, i.e., u^Xu = d. Then l>cp = J2^=i^jP^l- 
TZcp = {Rk} be the CP recovery map for $cp- Assume that 
p is in the code space, i.e., PpP = p. We now show that 
Ti-Cp[^L{p)] = A we have CP recovery. First, 



7^cp[$L(p)] 




N N 



5](i?fcG,F)p(PG],4). (22) 



Now, note that 

PGldP ■ 



u*^unPF}F,P = V <fcAy-u,;F 



dkSkiP- 



Then the polar decomposition yields 

GkP = C/fe(PGlGfcP)i/2 ^ y^f/.p. 
The recovery operation elements are given by 

Rk = ulPk, 



(23) 



(24) 



(25) 



where Pfe ?/fcPJ7fe. Therefore Pfe = GkPUl/^. This 
allows us to calculate the action of the fcth recovery operator 
on the hheiTorillll]: 



(26) 



RkGiP = uIpIGiP ^ul{UkPGll^)GiP 



'dkP. 



Therefore, 



7^Cp[*L(/5)] = 



N 

E 

■'=1 j-j 



2N 



X 



N 



P 

2pTr7t 



i=l 



(4jxA4P)p(Pv/d^(5, 



— pE •^^+«.^ 

i=l 



Next note that, using condition (l2ni(iii) and trace preservation 
by ^-l: 

PE'^E.P - 27|^P=^2Tr7tp = P^pfp,P = P 

i 

1 



Tr7''^ 



Hence, finally: 



ncp[Mp)] = p 



(28) 



(29) 



for any p in the codespace. ■ 

Note that $cp(p) need not be trace 
preserving: Tr[$cp(p)] = ^^AiE^iEjE^ + 
Y.ZiK^ E'i)p\, and while Y.f=iK^ E, = / if $l is 
trace preserving, we do not have conditions on X]i!=i^J^i 

andEf=ii?fi?.'- 

We define the class of "CP-recoverable linear noise maps" 
{*&cpr} as those <i>L for which CP recovery is always possi- 
ble. By Theorem m this includes all $l for which P can be 
found satisfying conditions (|2T]i(i)-(iii). However, these con- 
ditions are not necessary. 



B. Non-CP-recoverable linear noise maps 

We now define "non-CP-recoverable linear noise maps" 
{'i'nCPR} as those $L for which non-CP-recovery is always 
possible. Theorem |5] shows constructively that {$nCPR} in- 
cludes all linear noise maps (f>L for which P can be found 
satisfying only conditions (|2T]i(i) and (ii). Clearly, {<f>cp} C 

{*CPr} C {$nCPR} C {$l}- 

Theorem 5 Let <i>L = {Pi, E[}i be a linear noise map. Then 
every state p — PpP encoded using a QEC code defined by a 
projector P satisfying only Eqs. ( I27l )( i) and ( ii) can be recov- 
ered using a non-CP recovery map. 

Proof. Let Gk — Yl,i '^ikEi and GJ,, = Yl,i '^'ikE'i^ where the 
unitaries u and u' respectively diagonalize the Hermitian ma- 
trices a and a': d = u^au and d' = u'^^a'u' . Define a recov- 
ery map TZ — {Pfc, PJ,} (not necessarily CP) with operation 
elements 



Rk = UlPk, 



Rk ~ Uf} P^.. 



(30) 



^k-R^k are projection opera- 



(27) 



Here Pfc = UkPUl P[ 

tors, and Uk and U'^. arise from the polar decomposition of 
GkP and G'feP, i.e., GfeP = Uk{PGlGkP)'^^ and G'^P = 
UkiPG'^G'f.PY^'^ . The proof is entirely analogous to the 
proof of Theorem m except that we must keep track of both 
the primed and unprimed operators. Following through the 
same calculations we thus obtain RkGi^fp = ^/dk6kiy/p and 
R'k^iy/P ~ \f^k^ki\fp- Using this in the recovery map ap- 
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plied to the linear noise map, we find: 

nmPpP)] = Y^RkEiPpPE'i^R'^ 

kl 

= Y.Rk{Y.ul^G,)PpP{Y,u'uG'^)B!l 

kl j i 

= Fi^PpP^p, (31) 

where 

Fl = u*iju'i,^dkd'^5ko6M = Yl <ku'ik\l dkd'^ 

ijkl kl 

= Tr[w'd'tdMt] ^ Tr[u'Mtaa't] (32) 

is a "correction factor" for non-CP recovery of linear noise 
maps, which was 1 in the case of CP recovery, above. ■ 

Gathering the expressions derived in the last proof, we have 
the following explicit expressions for the left and right recov- 
ery operations: 

= UlPl = -^pY^u*,eI i?; = -j^pJ^uT.E'^. 

^ i V k i 

(33) 

This also shows that, in general, Rk need not equal i.e., 
the recovery map is linear but not necessarily CP. 

Note that standard QEC can also be interpreted as "error 
correction by inversion", in the following sense: when the 
noise map is CP and recovery is also CP, recovery is the in- 
verse of the noise map restricted to the code space (Theorem 
III. 3 in Ref. [7]). The same is true for our LQEC results above, 
which relax the restriction to CP noise maps. 

C. The physical case: Hermitian maps 

The general physical case is the case of Hermitian noise 
maps, to which any quantum dynamical process can be re- 
duced, as follows from Theorem [3] We can specialize Theo- 
rems|4]and|5]to this case. 

Corollary 1 Consider a Hermitian noise map $h (p) = 
^iLi^i^iP^i '^^d associate to it a CP map $cp(p) = 
X^il^i \ci\KipKl. Then any QEC code C and corresponding 
CP recovery map TZqp for <&CP <^fe also a QEC code and CP 
recovery map for $h- 

The important conclusion we can draw from Corollary[T]is 
that standard QEC techniques apply whether the noise map is 
CP or, as it will almost always be due to non-classical corre- 
lations, Hermitian. This is because Corollary [T] tells us that 
it is safe to replace all negative q coefficients by their abso- 
lute values, and thus replace the actual noise map by its CP 
counterpart. 

Proof. We have $h(p) = 'EZi^'^P^i with {E, = 
V^KijfLi and {E- = (y/ci)* Kij^^, whence we can apply 
the construction of Theorem ID Indeed, the "expanded" CP 

map becomes ^cvip) = ^Eti^.P^l + lEtiElpE'^ = 



Si=i \ci\KipKl , as claimed, and hence a QEC code and CP 
recovery for <I>cp is also a QEC code and CP recovery for <I>h. 
In particular, 7^cp[^h(/o)] = P- ■ 

Note that ^cp need not be trace preserving even in the Her- 
mitian map case: Tr[$cp(p)] = Tr^^i \c,\kIk,p], but if 
(^H is trace preserving then we only have 'Y^^=i^i-^l-^i ~ 
hence cannot conclude more about Tr[$cp(p)]- Also note 
that substitution of Ei — ^/clKi and E[ ~ {^fci)*Ki into the 

QEC conditions (|2B(i)-(iii) yields a-^ = 

and 7y = onj , i.e., unlike in the general linear maps 

case, the matrices a' and 7 in Eq. (l20l i are not independent 
from OL. In fact, as shown in Appendix|B]we can give a direct 
proof of Corollary [U which only invokes a single block of the 
A matrix. 



1. Example of CP recovery: Inverse bit-flip map 

Consider "diagonalizable maps", i.e., ^'d(p) = 
'Ylii CiKipKj, where Ci € C. The expanded CP map is 

$CP — J2i\'^i\-^iP-^i ■ Now consider as a specific in- 
stance an independent-errors inverse bit-flip map on three 
qubits: <I>ipf(p) = cop + ciJ2n=i ^npXn, where X„ is 
the Pauli ax matrix applied to qubit n, where co and ci 
are real, have opposite sign, and cq + 3ci = 1 (a Her- 
mitian map). Then l>cp = |co|p + \ci\Yfn=i ^^pXn, 
which is a non-trace preserving version of the well 
known independent-errors CP bit-flip map. The code 
is C = span{|OL) = |000),|1l> = |111)}, and 
P = |0l)(0l| + which satisfies Eq. dBTT ) 

with Fi ~ and ^2.3,4 = ■\/|ci|Xi_2,3. Then by 

Corollary [T] the same code (and corresponding CP recovery 
map) also corrects $ipf. The CP recovery map TZcp has 
operation elements Rq ^ P and = -^PXn}n=i'^ 

indeed, it is easily checked that TZcp[^ipf{P pP)] — PpP 
for any state p E C. 

2. Hermitian recovery maps 

Since Hermitian maps are the most general physical maps, 
it is natural to consider Hermitian recovery of Hermitian noise 
maps. We thus define "Hermitian recovery maps" {7?.h} as 
those Hermitian maps that correct a Hermitian noise map 
i.e., 7^H°*i'H(p) oc p. The following result presents a possible 
set of Hermitian recovery maps. 

Corollary 2 Consider a Hermitian noise map $h (p) = 
'Y^i^xCiKipK\ with error operators \Ki\ satisfying the re- 
lations PKJKjP — a.ijP. Any Hermitian map TZyi{p) = 
X]fe iT-kRkpRk with recovery operators {Rk} as in Eq. ( 1251 ) 
and {hk} € K corrects the noise map $h- 

The proof is given in Appendix ICl and employs a method 
similar to that of the proof of Theorem|5] 
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FIG. 1: The initial system-bath state is the generically non-VQD 
state psB{to)- The encoded system S — D + E consists of data 
qubits D and encoding ancillas E. We also include the recovery 
ancillas R, which are assumed to be completely isolated until they 
are brought into contact with S and B at a later time. Thus the 
full initial state is pss(io) CS> Pfl(io). The overall evolution is gov- 
erned by the unitary Usrb which acts on the system S, the bath B, 
and eventually the recovery ancillas R, and is denoted by the large 
grey box. The state of the data qubits is p|^) — TTE,B[psB{to)], 
a state which is as close as possible (by isolating the system) to the 
desired pure data state The state of each of the encoding ancil- 
las is p|o> = Tr£;',fl[psB(to)], a state which is as close as possible 
(again, by isolating the system) to the desired pure encoding ancilla 
state |0). Here Tr_E,s denotes a partial trace over all encoding an- 
cillas and the bath, Tr^'.s denotes a partial trace over all but one 
of the encoding ancillas, and the bath. Ideally, the encoding uni- 
tary Us is then applied to the encoded system. This is of course an 
idealization since in reality the encoding operation will not be a per- 
fect unitary; instead what is really applied is UsRB(ti,to), which 
is supposedly close to the ideal Us ® Ir ^S) Ib- Thus, after the 
encoding the total state is psRaiti) = UsRB{ti,to)[psB{to) ® 
Pfl(io)lt/gfls(ii, io) and the encoded system state is ps(ti) — 
T^iR,B[pSRB{ti)]- The system is then passed through the noise 
channel for the purpose of either computation or communication, 
i.e., pSRB{t2) = UsRB(t2,ti)psRB{ti)Uljig{t2,ti), whence 
Ps(fe) = TrH,s[psHs(f2)] = ^H[ps{ti)], where $h is a Her- 
mitian noise map since psRB(ti) is generically a non-VQD state 
due to the initial non-classical correlations between S and B. The 
goal of the error correction procedure is to recover the original en- 
coded system state from ps(i2), and to this end we introduce re- 
covery ancillas R at t2. Similarly to the encoding ancillas, these 
recovery ancillas are each in the state p|o) = Tts,r/ ,B[psRB{t2)], a 
state which is as close as possible to the desired pure recovery ancilla 
state |0). Next, ideally the recovery unitary Usr® Ib is applied. In 
reality what is applied is UsRB(T,t2), which is supposedly close 
to the ideal Usr ® Ib- Then the recovery ancillas are discarded 
and possibly recycled, leaving the encoded system in the final state 
Ps{T) = TrH,s[psijs(r)] = TZs[psit2)], which can be measured. 
Since psRB(t2) is generically not a VQD state (due to non-classical 
correlations between S and R, mediated by their mutual interaction 
with B), it is clear that the recovery map TZs is generically a non- 
CP Hermitian map. We recover the CP recovery map scenario if, 
for example, psRB{t2) ~ pSB{t2) ® pR{t2)- The assumption that 
this is not the case is consistent with the working premise of this pa- 
per and is equivalent in that regard to the assumption that the initial 
system-bath state is not of the form pss(to) = Ps(io) ® ps(to). 



3. How does non-CP, Hermitian recovery arise? 

In standard QEC theory the recovery map is considered CP. 
The reason for this is that the recovery ancillas are introduced 
after the action of the noise channel so that they enter in a ten- 
sor product state with the encoded qubits that underwent the 
noise channel. The recovery map is obtained in the standard 
setting by first applying a unitary over the encoded qubits plus 
recovery ancillas, then tracing out the recovery ancillas. This 
is manifestly a CP map over the encoded qubits. 

Since we know that the recovery map experienced by the 
encoded qubits is CP if and only if the initial state of the en- 
coded and recovery ancilla qubits has vanishing quantum dis- 
cord ifTsll . it is clear how a non-CP recovery map can be im- 
plemented: the recovery ancillas should have non-vanishing 
quantum discord with the encoded qubits. Since this will still 
be a QDP, the resulting recovery map will be Hermitian ac- 
cording to Theorem[3] 

Such a situation can come about in various ways. For ex- 
ample, a scenario which is particularly relevant for quantum 
computation and communication, is one where the environ- 
ment causes the recovery ancillas to become non-classically 
correlated with the encoded qubits before the recovery opera- 
tion can be applied. This is a reasonable scenario since, while 
the recovery ancillas are presumably kept pure and isolated 
from the environment for as long as possible, at some point 
they must be brought into contact with the encoded qubits, 
and at this point all qubits (encoded and recovery ancillas) are 
susceptible to correlations mediated by the environment. This 
is shown in Fig. [1] 



V. CONCLUSIONS 

This work aimed to fill two gaps: one in the theory of open 
quantum systems, and a resulting gap in the theory of quantum 
error correction. The first gap had to do with the type of maps 
that describe open systems given arbitrary initial states of the 
total system. In fact, it was not a priori clear that there should 
even be a linear map connecting the initial to the final open 
system state for arbitrary initial total system states. Build- 
ing upon the class of "special linear states" we introduced in 
II13I1 we showed here that in fact such a linear map description 
does always exist, and moreover, for quantum dynamics the 
map is always Hermitian. The map reduces to the completely 
positive type if and only if the initial total system state has 
vanishing quantum discord 1. 1 3.1 ; in all other cases it is Her- 
mitian but not CP. This result, we argued, impacts the theory 
of quantum error correction, where previously the assumption 
of CP maps was taken for granted. In the second part of this 
work we filled this gap in QEC theory, by developing a theory 
of Linear Quantum Error Correction (LQEC), which general- 
izes the CP-map-based standard theory of QEC. We showed 
that to every linear map $l is associated a CP map which, if 
coiTectable, also provides an encoding with coiTesponding CP 
recovery map for $l (Theorem|4|i. Moreover, it is possible to 
find a non-CP recovery for <f>L within a larger class of codes 
(Theorem |5]). From a physical standpoint this result is actu- 
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ally too general, since only Hermitian maps ever arise from 
quantum dynamics [to the extent that the standard quantum 
dynamical process (|2]i is valid]. Hence we specialized LQEC 
to the Hermitian maps case, and showed that in this case stan- 
dard QEC theory for CP maps akeady suffices, in the sense 
that it is legitimate to replace a given Hermitian noise map by 
a corresponding CP map obtained simply by taking the abso- 
lute values of all the Hermitian map coefficients. Any QEC 
code which corrects this CP map will also correct the original 
Hermitian map (Corollary [1]). Nevertheless, there is room for 
a genuine generalization when one considers Hermitian maps, 
since it is also possible to perform QEC using Hermitian re- 
covery maps (Corollary |2]i. We argued that, in fact, recov- 
ery maps will generically be non-CP Hermitian maps, since 
recovery ancillas that are introduced into a quantum circuit 
prior to the recovery step will become non-classically corre- 
lated with the environment and consequently with the rest of 
the system. 

An interesting open question for future studies is whether 
the results presented here have an impact on the threshold for 
fault tolerant quantum error correction. For example, note that 
while CP recovery perfectly returns the encoded state [Eqs. 
( |29] l and (IB8I 1I. non-CP recovery only does so up to a propor- 
tionality factor which depends on the details of the noise and 
recovery maps [Fl in Eq. ( [32] i and Fh in Eq. ( |C2| )1. This pro- 
portionality factor - assuming non-CP recovery is applied - 
may differ for different terms in the fault path decomposition 
112811 . an effect which may propagate into the value of the fault 
tolerance threshold. This requires careful analysis, which is 
beyond the scope of this paper. 
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APPENDIX A: PROOF OF THEOREM 1 

We use a method similar to Choi's proof for a CP map 
representation [34], recently clearly reviewed in Ref. Il35ll . 
The main difference between the proofs in Refs. ll34ll35ll and 
our proof is that in the previous proofs positivity allowed for 
the use of standard diagonalization, whereas in the absence of 
positivity we use the singular value decomposition I3Q1 . 
Proof. Eq. (|4| immediately implies that $l is a linear map. 
For the other direction, let AI = X]"j=i K)OI ® l*)OI = 
n\(j)) (01, where \i) is a column vector with 1 at position i and 
O's elsewhere, and \(f)) = rT^I"^ |i) ® \i) is a maximally 
entangled state over H (E) H, where H is the Hilbert space 
spanned by M is also an n x n array of n x n ma- 

trices, whose (i, j)th block is Construct two equivalent 



expressions for (T® $l) [Af], where T is the {nxn)x {nx n) 
identity matrix, (i) {T (g) is an n x n array of m x m 

matrices, whose {i,j)th block is $L[|i)0|]- (ii) Consider 
a singular value decomposition: (X (g) $)[M] — UDV — 
Y.a KU\a){a\V = J2a >^a\ua){va\- Here U and V are uni- 
tary, D = diagdAo,}) is diagonal and Aq > are the singular 

values of (J (g) Divide the column (row) vector \ua) 

{{va I) into n segments each of length in and define an to x 71 
{n X m) matrix Ea {E'^) whose ith column (row) is the ith 
segment; then Ea\i) {{i\E'^) is the ith segment of \ua) {{va\)- 
Therefore the (i,j)th block of |uQ)(fQ| becomes Faji) 

Equating the two expressions in (i) and (ii) for the 

(i,j)th block of (X (g $l)[M], we find $L[|«)(j|] = 
'Yl,a'^aEa\i){j\E'^. Sincc Ac > we can redefine Ea as 
\/\^Ea and E'^ as \/\^E'^, which we do from now on. Fi- 
nally, the linearity assumption on <I>l, together with the fact 
that the set spans 9Jl„, implies Eq. 

Next let us prove Eq. ^ for Hermitian maps. For an old 
proof that uses very different techniques see Ref. iIstIi . Eq. Q 
immediately implies that $h is a Hermitian map. For the other 
direction, associate a matrix L,^^ with the Hermitian map 
p' — <f>H(/o) p'mii — ^nvPnv (summation ovcr repeated 
indices is implied). Hermiticity of p and its image p' implies 

„' — „'* Tmti* * rmfj.* j „ Tmfi* _ j fim 

Putn — Prnfi ^ ^nv Pnu ~ ^nu Pi"" ^nu ~ ^un 

[38]. We can use this property of L,^,^ to show that if <I>h is a 
Hermitian map, then I ® <i>H is Hermiticity preserving. Con- 
siderX = Mll\k){£,\®\n){v\. Then M' = {I®<^yi)[M] = 
Mll\k){^\ ® $H(|n)(H) - Mr^nfc)(^l ® LZ^HH As- 
sume that = Al^". This property holds for M = 

M = where |0) = dim{n)~^^'^ J2i K) ® N) is a maxi- 

mally entangled state over H (g 7^ (A^^™ = 1). Then M'^ = 

M^nOikl ® LZ^*W){n\ = Mir\0{k\ ® k>(n| = 
Ai'. Therefore (I (g) <I>h)[|(?!')(0|] is Hermitian, and in par- 
ticular invertible. It follows that the SVD used in the proof 
of Theorem [T] can be replaced by standard diagonalization 
(U = V^). In this case the left and right singular vectors 
\ua) — {va\^ are the eigenvectors of (X ® $h)[|0)(0|] and 
Ca = Aq are its eigenvalues. Then Ea — E'a in Eq. (|4|i and 
Ca e M. ■ 

We note that by splitting the spectrum of {T ® $h) \\4>) ('/'I] 
into positive and negative eigenvalues, {c+ > 0} and {c^ < 
0}, we have as an immediate corollary a fact that was also 
noted in llloll : Any Hermitian map can be represented as 
the difference of two CP maps: $(p) = '^a^a^aP^a^ ~ 

EaKlKpE-^- 



APPENDIX B: DIRECT PROOF OF COROLLARY 1 

Proof. The operation elements of $cp are {Fi = 
^\K,}fL„ whence I>cp(/9) = EZlF^pFl The standard 
quantum error conditions (fT9b for <&cp is a set of conditions 
in terms of the Fi : 

PFIf,P = I3,,P, €{!,..., TV}. (Bl) 
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The existence of a projector P which satisfies Eq. ( IBll ) is 
equivalent to the existence of a QEC code for <i>cp ■ Assum- 
ing that a code C has been found (i.e., PC = C) for $cp, we 
use this as a code for $h and show that the corresponding CP 
recovery map 7?.cp is also a recovery map for <i>H- Indeed, let 
Gj = X^ili "jj^i be new operation elements for $cPi i-S-^ 
"^cp — SjLi^jP^]' where u is the unitary matrix that di- 
agonalizes the Hermitian matrix /3 = [f^ij], i.e., u^(3u — d. 
Let 7?,cp = {Rk} be the CP recovery map for $cp- Assume 
that p is in the code space, i.e., PpP — p. We now show that 
T^CP [*I'h(p)] = P, i-e., we have CP recovery. First, 

7^cp[<^>H(p)] - E^fc 1^^^/^^^) 

N N 

= Eui E 

x^(i?,.G,P)p(p44). (B2) 

k 

Now, note that, using Eq. ( IBll l: 

(B3) 



Then the polar decomposition yields GkP — 
UkiPGlCkPy^ = VchUkP. The recovery operation 
elements are given by 



Rk = UlPk-, Pk = UkPUl 



(B4) 



Therefore Pk — GkPUl/^/dk. This allows us to calculate 
the action of the fcth recovery operator on the Ith error: 



RkGiP = uIpIGiP = uliUkPGl/ ^)GiP 



'dkP. 



(B5) 



Therefore, 

7^cp[*H(p)] 



N N 

c 



E^ E 

■'=1 jj'=i 

^ E i^^kjVdkP^ P {P\fdk6uj^ 



k 
N 

N 



= (E ih^^^P' 



(B6) 



Next note that, using condition ( IBll ) and trace preservation by 
$h: 



N 



^=l 



^E O^^^'^ ^ Pj2c,KjK,P = P 



TV 



Hence, finally: 



Z— 1 ' ' 



7ecp[$L(p)] = p 



(B7) 



(B8) 



for any p in the codespace. 



APPENDIX C: PROOF OF COROLLARY 2 

Proof. Let {P,; = yJ\ci\Ki\f^^\ we simply use the identities 
given in the proof of the previous theorem - specifically Eq. 
- to calculate TZ-a o <I'h(p) 



N ^ N 



N 



= PpP^--!-^hkU*kU^kdk 
= F^iPpP^p, 



(CI) 



where 



(C2) 



where h = diag({ft,fc}), and Ph is a "correction factor" for 
Hermitian recovery of Hermitian noise maps, which was 1 in 
the case of CP recovery, above. ■ 



[1] E.B. Davies, Quantum Theory of Open Systems (Academic [2] R. Alicki and K. Lendi, Quantum Dynamical Semigroups and 
Press, London, 1976). Applications, No. 286 in Lecture Notes in Physics (Springer- 



11 



Verlag, Berlin, 1987). 

[3] H.-P. Breuer and F. Petruccione, The Theory of Open Quantum 
Systems (Oxford University Press, Oxford, 2002). 

[4] M.A. Nielsen and l.L. Chuang, Quantum Computation and 
Quantum Information (Cambridge University Press, Cam- 
bridge, UK, 2000). 

[5] F. Gaitan, Quantum Error Correction and Fault Tolerant Quan- 
tum Computing (CRC Press, Boca Raton, 2008). 

[6] K. Kraus, States, Effects and Operations, Fundamental Notions 
of Quantum Theory (Academic, Berlin, 1983). 

[7] E. Knill and R. Laflamme, Phys. Rev. A 55, 900 (1997). 

[8] R. Alicki, Phys. Rev. Lett. 75, 3020 (1995), Comment on "Re- 
duced Dynamics Need Not Be Completely Positive"; Reply by 
P. Pechukas, ibid., p. 3021. 

[9] P Stelmachovic and V. Buzek, Phys. Rev. A 64, 062106 (2001). 
[10] T.F Jordan, A. Shaji, and E.C.G. Sudarshan, Phys. Rev. A 70, 
052110(2004). 

[11] H. A. Carteret, D. R. Temo, and K. Zyczkowski, Phys. Rev. A 

77, 042113 (2008). 
[12] C.A. Rodriguez, K. Modi, A.-M. Kuah, E.C.G. Sudarshan, A. 

Shaji, J. Phys. A 41, 205301 (2008). 
[13] A. Shabani and D. A. Lidar, eprint arXiv.org:0808.0175 (in 

press, Phys. Rev. Lett.). 
[14] H. GUivier and W.H. Zurek, Phys. Rev. Lett. 88, 017901 (2001). 
[15] R. Alicki, D. A. Lidar, and R Zanardi, Phys. Rev. A 73, 052311 

(2006). 

[16] D. M. Tong, L. C. Kwek, C. H. Oh, J.-L. Chen, and L. Ma, 

Phys. Rev. A 69, 054102 (2004). 
[17] W. F Stinespring, Proc. Am. Math. Soc. 6, 211 (1955). 
[18] L. Jakobczyk and M. Siennicki, Phys. Lett. A 286, 383 (2001). 
[19] G. Kimura, Phys. Lett. A 314, 339 (2003). 
[20] M.S. Byrd and N. Khaneja, Phys. Rev. A 68, 062322 (2003). 



[21] G. Kimura and A. Kossakowski, Open Sys. and Information 
Dyn. 12, 207 (2005). 

[22] P.W. Shor, in Proceedings of the 37th Symposium on Founda- 
tions of Computing (IEEE Computer Society Press, Los Alami- 
tos, CA, 1996), p. 56. 

[23] D. Aharonov and M. Ben-Or, in Proceedings of 29th Annual 
ACM Symposium on Theory of Computing (STOC) (ACM, New 
York, NY, 1997), p. 46. 

[24] E. Knill, R. Laflamme and W. Zurek, Science 279, 342 (1998). 

[25] A.M. Steane, Phys. Rev. A 68, 042322 (2003). 

[26] E. Knill, Nature 434, 39 (2005). 

[27] B.W. Reichardt, in Automata, Languages and Programming, 
Vol. 4051 of Lecture Notes in Computer Science, edited by M. 
Bugliesi et al. (Springer, Berlin, 2006), Chap. 6, p. 50. 

[28] D. Aharonov and M. Ben-Or, SIAM J. Comput. 38, 1207 

(2008). 

[29] P Aliferis and J. Preskill, Phys. Rev. A 79, 012332 (2009). 

[30] B.M. Terhal, G. Burkard, Phys. Rev. A 71, 012336 (2005). 

[31] P. Aliferis, D. Gottesman, and J. Preskill, Quantum Inf. Com- 
put. 6, 97 (2006). 

[32] D. Aharonov, A. Kitaev, J. Preskill, Phys. Rev. Lett. 96, 050504 
(2006). 

[33] M.A. Nielsen, CM. Caves, B. Schumacher and H. Bamum, 
Proc. Roy. Soc. London Sen A 454, 277 (1998). 

[34] M.D. Choi, Linear Algebr. Appl. 10, 285 (1975). 

[35] D.W. Leung, J. Math. Phys. 44, 528 (2003). 

[36] R.A. Horn and C.R. Johnson, Matrix Analysis (Cambridge Uni- 
versity Press, Cambride, UK, 1999). 

[37] R.D. Hill, Linear Algebra and its Applications 6, 257 (1973). 

[38] K. Zyczkowsld and 1. Bengtsson, Open Sys. & Information 
Dyn. 11, 342 (2004). 



